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Abstract 

A wind  tunnel  interference  assessment  concept 
which  presents  a rational  predictive  means  of  wall 
interference  analysis  is  evaluated.  The  procedure 
consists  of  employing  as  an  outer  boundary  condi- 
tion an  experimentally-measured  pressure  distribu- 
tion along  a convenient  control  surface  located 
inward  from  the  actual  tunnel  walls.  Attention  has 
been  focused  on  axisymmetric  flows  in  the  transonic 
regime  where  tunnel  interference  is  high  and  where 
the  experimentally-measured  conditions  on  the 
control  surface  are  of  mixed  subsonic/supersonic 
type.  Based  on  the  transonic  small-disturbance 
equation,  results  for  surface  and  near  flow  field 
pressure  distributions  are  presented  for  a variety 
of  different  slender  body  shapes.  These  calcula- 
tions indicate  both  the  accuracy  of  the  procedure 
as  well  as  its  ease  of  implementation.  The  pro- 
cedure relates  directly  to  the  correctable- 
interference  wind  tunnel  concept  recently  suggested. 


Introduction 

The  assessment  of  transonic  wind  tunnel  inter- 
ference remains  an  unsatisfactorily  resolved 
problem  after  almost  three  decades  use  of  the  venti- 
lated tunnel  wall  concept.  The  classical  homo- 
geneous wall  boundary  conditions^  for  modeling 
perforated  or  porous  and  slotted  walls  have  proven 
to  be  of  considerable  value  for  providing  insight 

into  the  overall  effects  of  ventilated  tunnel  walls 
2 3 

at  transonic  speeds  ' , but  as  presently  consti- 
tuted, they  are  seriously  deficient  for  a quanti- 
tative study  of  wall  interference.  Although  the 
wall  characteristics  are  known  to  depend  on 
geometrical  (test  article,  tunnel  ventilation 
geometry)  as  well  as  flow  parameters  (M^,,  Re, 
local  blowing  pressure) , the  exact  nature  of  the 
wall  flow  is  not  adequately  known.  It  has  been 
widely  assumed  that  it  is  basically  viscous  in 
character4,  although  some  recent  slotted  wall 
results  by  Berndt^  indicate  that  under  certain  con- 
ditions the  wall  flow  for  slotted  walls  may  not  be 
completely  dominated  by  viscosity.  For  porous 
walls,  it  has  been  found6,7,  that  the  wall  boundary 
layer  can  create  a nonlinear  relationship  between 
velocity  normal  to  the  wall  and  the  difference  in 
pressure  across  the  wall.  Kacprzynski8  has 
demonstrated  the  feasibility  of  two-dimensional 
potential  flow  calculations  using  an  assumed  non- 
linear wall  characteristic.  The  present  state  of 
development,  however,  of  ventilated  wall  boundary 
conditions  is  such  that  they  are  inadequate  as  a 
predictive  method;  namely,  they  cannot  provide  a 


quantitative  a priori  assessment  of  tunnel  inter- 
ference at  transonic  speeds. 

The  adaptive  wall  concept  is  an  attempt  to 
remedy  this  situation  by  combining  an  interference 
assessment  capability  together  with  a tunnel  wall 
modification  capability  to  actually  eliminate  the 
interference.  Notwithstanding  the  problems  associ- 
ated with  the  implementation  of  adaptive  wall 
technology,  the  fact  remains  that  many  of  the 
currently  operating  transonic  tunnels  possess  only 
a limited  degree  of  wall  control  or  none  at  all. 
Consequently,  the  concept  of  a correctable-inter- 
ference transonic  tunnel  recently  proposed9,  and 
toward  which  the  procedures  described  in  this  paper 
are  directly  related,  appears  to  have  strong  merit. 
The  heart  of  the  correctable-interference  concept 
is  that  the  distribution  of  a single  experimentally- 
measured  flow  quantity  (e.g.,  pressure  or  flow 
angle) , obtained  along  a convenient  control  surface 
located  sufficiently  inward  from  the  actual  tunnel 
walls  so  as  to  be  removed  from  local  wall  dis- 
turbances, is  employed  as  an  outer  boundary 
condition.  This  is  then  used  together  with  a 
calculative  procedure  to  determine  the  potential 
flow  about  the  model  interior  to  the  control 
surface.  Comparison  of  calculated  and  measured 
flow  quantities  at  the  model  surface  serve  to 
assess  the  influence  of  nonpotential  effects 
(viscosity,  rotation) ; while  comparison  with  a 
calculated  free  air  flow  can  provide  a quantitative 
rational  assessment  of  tunnel  interference  effects. 

In  this  paper  we  formulate  the  assessment 
procedure  based  on  a measured  pressure  distribution 
boundary  condition,  describe  its  application  to 
axisymmetric  transonic  flows  past  a variety  of 
slender  bodies,  and  demonstrate  its  effectiveness 
by  extensive  comparisons  with  data.  The  theoretical 
predictions  employ  finite-difference  SLOR  solutions 
of  the  axisymmetric  transonic  small -disturbance 
potential  equation.  Results  are  presented  for 
pressure  distributions  on  the  surface  and  in  the 
near  flow  field  of  these  various  slender  bodies  at 
free  stream  Mach  numbers  throughout  the  transonic 
regime.  In  order  to  provide  a severe  test  of  the 
procedure,  the  particular  geometries  and  flow  con- 
ditions were  purposely  selected  to  be  in  a range  of 
free  stream  Mach  numbers  near  one  where  the  tunnel 
interference  is  most  pronounced,  and  where  the 
experimentally-measured  data  on  the  control  surface 
are  of  mixed  subsonic/ super sonic  character. 
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JL  J and  < > signify  the  difference  and  the  mean, 
respectively,  of  the  enclosed  quantity  on  the  two 
sides  of  the  shock  surface. 


Analysis 


Basic  Equations 

The  coordinate  system  employed  in  the  analysis 
is  a body-fixed  cylindrical  system  with  origin  at 
the  nose  of  the  body  as  illustrated  in  figure  1. 


VENTILATED  TUNNEL  WALL  (r  - H/t) 


CONTROL  SURFACE  (r  - r^) 


Fig.  1 Illustration  of  model  and  tunnel  geometry 

The  flow  is  assumed  to  be  inviscid  and  steady,  and 
the  body  shapes  sufficiently  slender  and  smooth 
that  the  resulting  flow  field  is  irrotational  and 
adequately  treated  by  small-disturbance  theory. 
Accordingly,  a disturbance  velocity  potential  $ can 
be  defined  by 

$(x,r)  = [x  + $ (x,r)  ] (1) 

where  * is  the  total  velocity  potential,  U„  repre- 
sents the  uniform  oncoming  velocity,  t is  the  com- 
plete body  length,  and  the  coordinates  (x,r)  have 
been  nondimensionalized  by  t.  The  governing 
partial  differential  equation  for  $ is  the  usual 
axisymmetric  transonic  small-disturbance  equation 

(1  - M2  - l/2M^(y  + l)$x  + p (r*r)r=  0 (2) 

The  pressure  coefficient  in  the  vicinity  of  the 
slender  body  is  given  by 

Cp(x,r)  = -24>x(x,r)  - <)>2(x,r)  (3) 

With  regard  to  the  flow  domain  indicated  in  figure 
1,  the  boundary  conditions  to  be  imposed  on  the 
solution  consist  of  ( i)  outer  flow  conditions  on 
the  upstream,  downstream  and  lateral  boundaries 
which  are  appropriate  to  the  behavior  of  the  body 
in  a free-air  or  wind  tunnel  environment,  (ii)  the 
body  surface  condition  of  no  normal  flow,  and  (iii) 
shock  wave  conditions  to  be  applied  at  any  shock 
surface  appearing  in  the  flow  such  that  the  poten- 
tial is  continuous  through  the  shock  and  the 
velocity  components  satisfy  the  small-disturbance 
approximation  to  the  Rankine-Hugoniot  conditions 
at  the  location  of  the  shock.  The  appropriate 
outer  flow  boundary  conditions  are  discussed  in  the 
following  section.  The  requirements  on  the  body 
and  shock  surface  lead  to  the  following  conditions 
on  t 


lim  (r$r) 
r-*o 


S'  (x) 
2n 


(4) 


I ♦ a, hock  ■ 0;  [rU-M2-M2(y+l)<V}x 
+ ^ 3>]  shock  " 0 


where  S (x)  is  the  body  cross  sectional  area  distri- 
bution nondimensionalized  by  l 2,  and  the  symbols 


Outer  Boundary  Conditions 

The  outer  boundary  condition  to  be  applied 
along  the  cylindrical  control  surface  indicated 
in  figure  1 is  determined  from  static  pressure 
measurements.  The  actual  condition  imposed  on  C is 
of  Dirichlet  type,  and  is  given  by 


♦ <x,rg)  = -1/2  I Cpm(5,rs)dt  + 4>  (x^r^  (6) 

where  r denotes  the  radial  location  of  the  control 
surface,  x^  the  position  of  the  upstream  boundary, 
and  Cpm  the  measured  pressure  coefficient  along 
rs.  The  quantity  if(x^,rs),  which  is  proportional 
(to  within  a constant)  to  the  average  flow  inclina- 
tion at  the  upstream  boundary,  has  for  convenience 
and  without  loss  of  generality  been  set  to  zero. 
Along  the  upstream  and  downstream  boundaries,  for 
subsonic  oncoming  flow  we  have  employed  the  condi- 
tions 

< 1:  $(x1(r)  = 0;  4>(xN,r)  * ♦(xN,rs)  (7) 

while  for  sonic  and  supersonic  oncoming  flow,  we 
have  used  the  uniform  flow  conditions 

M«,  >.  1:  (Mx^r)  = ♦x*xl»r)  = 0;  <(>x(xN,r)  = 0 


While,  strictly  speaking,  for  Moo  > 1,  no  condition 
is  required  at  the  downstream  boundary  if  as  we 
assume  it  is  sufficiently  far  removed  so  that  the 
outflow  is  entirely  supersonic.  For  actual 
numerical  application,  however,  some  condition  is 
needed*  to  treat  subsonic  outflow  which  may  develop 
during  the  course  of  the  relaxation  solution 
process  prior  to  convergence.  From  numerical 
experimentation,  the  condition  $x  = 0 given  in 
equation  (8)  has  been  found  to  be  satisfactory. 
These  conditions  are  summarized  in  figure  2. 
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Fig.  2 Summary  of  boundary  conditions  for  free 
air  and  wind  tunnel  environment 


For  the  corresponding  comparative  free  air 
calculations,  we  have  employed  the  following  condi- 
tions. For  Moo  < 1 flows,  the  asymptotic  far-field 
solution  4ff  given  by  an  axisymmetric  source  and 
axisymmetric  doublet*0  is  employed  on  all  the  outer 
boundaries,  i.e. 


2 


♦jj*  <x,r) 


if"  + 1/4  l’S 


S(C)d5  + 


M^ty  + 1)  J J *2<t,r)rdrdU  • ^ (9) 

-00  o r 

where  r - (x2  + (1  - M£)r2)1/2,  x > x - 1/2,  xq  is 
the  location  of  the  body  base,  and  SQ  is  either  the 
base  area  of  the  body  or  that  of  a sting  attached 
to  the  body  base  nondimensionalized  by  l2.  On  the 
lateral  boundary,  for  M„  = 1 flows,  we  employ  the 
asymptotic  solution** 


(x,r) 


f(x/r4/7) 


while  for  > 1,  the  following  relation*  is  used 


1 (x,r)  - f (x  - /m2-1  r- 


M|(y+1)  • r*x 

/m2  - 1 


On  the  upstream  and  downstream  boundaries,  for 
>.  1,  the  conditions  given  in  equation  (8)  are 
employed.  These  conditions  are  also  sunmarized  in 
figure  2. 

Computational  Procedure 

The  method  employed  for  the  solution  of 
equation  (2)  is  a finite-difference  successive  line 
over- relaxation  (SLOR)  procedure  using  Murman-Cole 
type-dependent  difference  operators1”3.  To  realize 
the  calculation,  we  have  employed  the  following 
fully- conservative  form 

[n(K*  - 1/2?  2)J  + In?  1 = o (12) 

x xx  n n 


*(x,n)  - ♦(x,r)/T2;  K - (1-M2) /t2M2 (y+1) ; 

n - T M^/y+1  r (13) 

and  t denotes  the  thickness  ratio  of  the  body.  The 
finite  difference  form  of  the  equation  actually 
solved  is  that  suggested  by  Jameson**  in  terms  of  a 
correction  potentia}.  Additionally,  a pseudo-time 
term  of  the  form  -e—  4xt  was  added  to  enhance 
stability  and  speed  convergence. 


In  order  to  examine  the  feasibility  of  the 
procedure  described  above  for  assessing  tunnel 
interference,  as  well  as  to  examine  the  stability 
and  convergence  characteristics  of  the  computation- 
al method  subject  to  the  experimentally- imposed 
Dirichlet  condition  of  equation  (6) , we  have 
extensively  tested  the  method  using  data  obtained 
in  a conventional  transonic  tunnel  on  five  dif- 
ferent slender  body  shapes  and  at  Mach  numbers  at 
and  near  one.  The  particular  geometries  and  flow 
conditions  were  purposely  selected  to  study  situa- 
tions where  tunnel  interference  is  most  pronounced, 
and  insofar  as  possible  where  the  experimentally- 
measured  data  on  the  control  surface  are  of  mixed 
subsonic/supersonic  character. 


The  body  shapes  examined  are  (i)  parabolic- 
arc*®,  (ii)  power-law  body  with  SmaLX  at  x * 0.3*®, 
(iii)  power-law  body  with  S,^,,  at  x = 0.7*®,  (iv) 
parabolic-arc  with  bumpy  midsection*®,  and  (v)  para- 
bolic-arc with  indented  midsection*®.  Geometric 
details  of  these  shapes  are  provided  in  figure  3. 
Data  for  all  of  these  bodies  were  obtained  in  the 
Ames  14-foot  transonic  wind  tunnel. 
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Fig.  3 Geometric  details  of  various  wind  tunnel 

models  considered  with  dimensions  in  inches 

Smooth  Bodies 

In  figure  4 we  have  displayed  comparisons  of 
the  theoretical  results  with  data*®  for  a para- 
bolic-arc body  of  revolution  having  a diameter  to 
complete  body  length  D/£*=12.  Results  are  shown  at 
free-stream  Mach  numbers  M„  - 0.975,  1.00,  1.025, 
and  1.10.  For  the  particular  body  tested  the 
tunnel  half-height  to  complete  body  length  ratio 
was  H/£»7/6.  The  theoretical  results  indicated  by 

the  solid  lines  ( ) represent  SI/3R  solutions  of 

equation  (12)  employing  the  measured  pressure  dis- 
tribution outer  boundary  condition  equation  (6) . 
Pressure  survey  data  obtained  at  the  radial  loca- 
tion r/I>>4,  which  represents  the  outermost  survey 
station  at  which  data  were  taken,  were  used.  The 
results  indicated  by  the  dashed  lines  ( ) repre- 

sent the  corresponding  theoretical  results  for 
free  air  conditions,  equations  (9-11). 

The  theoretical  body  surface  pressure  coef- 
ficients were  determined  by  extrapolating  the  SLOR 
results  obtained  along  the  first  radial  grid  line 
r^  down  to  the  body  surface  according  to  the 
slender  body  result 

Cp ( x , R)  - - 2*x(x.r1)  + 025l«n(i)  ♦ (§) 


3 


where  R is  the  equivalent  body  radius.  The  flow 
field  pressures  were  calculated  by  bi-linear  inter- 
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Fig.  4 Experimental*®  ( q ) and  theoretical  (■ 
pressure  distributions  for  a parabolic-arc 
body  with  D/£-l/12  at  several  transonic 

Mach  numbers;  Cp^r/D-4)  input, 

free  air 


SLOR  solution  process,  relaxation  factors 
1,7  < u <1,9  were  used  and  the  coefficient  of 
the  pseudo-time  term  was  set  to  e - 0.5.  No 
stability  problems  were  encountered  and  convergence 
in  all  of  the  imposed  pressure  distribution  cases 
was  very  rapid,  with  a |A*|lnax  < 10"®  criteria 
reached  within  75  iteration  sweeps.  Pressure 
changes  on  the  body  usually  become  less  than  10“^ 
before  40  iterations.  The  corresponding  results 
for  the  free  air  calculations  generally  required 
about  25%  more  iteration  cycles  to  reach  the  same 
convergence  level . 

With  regard  to  the  comparisons  in  figure  4,  we 
note  that  the  agreement  between  the  data  and  the 
imposed  pressure-condition  theoretical  results  is 
excellent  for  both  surface  and  flow  field  pressures. 
The  location  and  strengths  of  shocks  are  accurately 
predicted  throughout  this  flow-sensitive  Mach 
number  range.  Viscous  effects  on  these  flows, 
which  have  Reynolds  numbers  based  on  body  length 
of  approximately  24  million,  appear  to  be  quite 
small  and  confined  to  the  flow  imnediately  behind 
the  body  base  where  the  step-down  sting  results 
in  a separated  flow  situation.  The  corresponding 
free  air  results  are  interesting  in  that  they 
indicate  essentially  no  tunnel  interference  at 
M«,  - 0.975,  a substantial  amount  at  M*,  - 1.00  and 
1.025,  and  then  very  little  again  at  M„  = 1.10. 

This  serves  to  point  out  the  narrowness  of  the 
critical  range  for  this  class  of  smooth  slender 
bodies. 

Figure  5 presents  analogous  results  for  a 
slender  body  having  the  same  maximum  diameter  to 
length  ratio  (D/£«l/12)  as  the  parabolic-arc  body 
of  figure  4,  but  with  the  location  of  maximum 
diameter  at  30%  of  the  body  length.  The  ordinates 
of  this  body  are  given  by 

R - A [(1  - x)  - (1  - x)nJ  (15) 

where  n «=  6.03  and  A is  related  to  D/l  and  n by 


polation  for  $ through  the  flow  field  grid,  and 
then  employing  equation  (3)  to  calculate  the  pres- 
sure coefficient.  The  (x,r)  mesh  density  used  for 
the  imposed  pressure  distribution  calculations  was 
(80  < 24)  with  40  equally-spaced  points  on  the 
body.  The  r-grid  as  well  as  the  x-grid  ahead  and 
behind  the  body  were  expanded  geometrically  using 
a grid  ratio  of  1.2:1.  The  x-mesh  extended  2 body 
lengths  ahead  and  2 body  lengths  behind  the  body. 
The  location  of  the  first  r-grid  line  was  at 
r = D/2 l and  the  control  surface  was  at  rs  - 1/3. 
For  the  free  air  calculations,  an  identical  x-mesh 
was  used  and  a 40  point  r-grid  was  employed  that 
was  identical  to  the  previous  r-mesh  out  to  r - 1/3 
but  which  was  continued  laterally  out  to  r * 5.  In 
nunerically  implementing  the  imposed  pressure  dis- 
tribution boundary  condition  equation  (6) , it  was 
necessary  to  extend  the  data15'16  to  the  upstream 
and  downstream  boundaries,  since  the  experimental 
pressures  were  only  obtained  for  -0.15  < x/l  < 1.0. 
In  addition,  numerical  experimentation  indicated 
that  some  degree  of  data  point  smoothing  was 
desirable.  Because  the  data  sometimes  contained 
one  or  more  shock  waves,  however,  and  the  faithful 
reproduction  of  those  shock  profiles  was  considered 
essential,  the  method  selected  was  a least  squares 
smoothing  spline  fit  with  individual  point  weigh- 
ing. As  can  be  seen  in  figure  4 as  well  as  in  the 
results  to  follow,  this  produced  extremely  satis- 
factory results  in  all  cases  considered.  In  the 


A * jyv'ta-D/j  (n— 1) ] (D/£)  (16) 

As  in  the  case  of  the  parabolic-arc  body,  compari- 
sons of  the  Cj^-input  theoretical  results  with  data 
again  display  excellent  agreement.  The  correspond- 
ing free  air  comparisons  indicate  that  the  tunnel 
interference  effects  are  more  pronounced  for  this 
shape  than  for  the  parabolic-arc.  This  is  to  be 
anticipated  and  is  in  accord  with  the  observations 
noted  in  reference  17  that,  for  a given  thickness 
ratio  and  length,  wave  reflection  interference  on 
the  afterbody  region  increases  as  the  point  of 
maximum  thickness  moves  forward. 

We  note  for  this  body  that  the  » 0.975  re- 
sults display  only  slight  interference,  while  the 
Mo,  - 1.00  comparisons  indicate  strong  tunnel  inter- 
ference which  has  resulted  in  a shock  movement  on 
the  body  surface  of  approximately  15%  of  the  body 
length.  For  the  M„  • 1.025  results,  the  measured 
boundary  condition  was  imposed  at  r/D  « 3 rather 
than  4,  since  r/D  - 4 data  were  not  obtained  for 
this  flow.  We  observe  that  the  comparison  of  the 
predicted  result  with  data  is  quite  good,  in  par- 
ticular, even  for  the  flow  field  pressures  at  x 
locations  behind  the  body  base.  Only  for  the 
Kn  » 1.10  case  does  a slight  discrepancy  arise,  and 
that  only  for  the  flow  field  results  at  locations 


4 


tunnel  interference  for  this  shape  is  very  small 
and  restricted  essentially  to  the  Mb  - 1.00 
results  at  axial  locations  behind  x » 0.90. 


SURFACE 


Fig.  6 Experimental'1  (§)  and  theoretical  (• 
pressure  distributions  for  a body  with  maxi 
mun  thickness  at  x/£= 0.7  and  D/i=l/12  at 
several  transonic  Mach  numbers) 

Cp^lr/DF-4)  input, 

free  air 


Fig.  5 Experimental15  (§>  and  theoretical  ( ) 

pressure  distributions  for  a body  With  max- 
imum thickness  at  x/Z=0.3  and  D/Z=l/12 
at  several  transonic  Mach  numbers: 

Cpjjfr/D-4)  input, free  air 


Bumpy  and  Indented  Bodies 


behind  the  body  base.  This  may  be  caused  by  the 
interaction  of  Em  oblique  shock  emanating  from  the 
body  base  and  interacting  with  the  separated  base 
flow.  A more  modest  indication  of  this  can  be 
observed  in  the  corresponding  results  at  M„  “ 1.10 
for  the  parabolic-arc  body  shown  in  figure  4. 


In  order  to  provide  an  indication  of  the 
ability  of  the  method  to  handle  situations  where 
multiple  shocks  or  several  high  gradient  regions 
exist  on  the  control  surface,  we  have  examined  the 
bumpy  and  indented  bodies16  which  were  tested 
transonically  in  the  Ames  14-foot  wind  tunnel. 

These  bodies  were  quite  long  (£  « 112  inches  - see 
figure  3) , and  consisted  of  basic  parabolic-arc 
bodies  having  maximum  diameter  to  length  ratios 
D/Z  *»  14  but  with  a sinusoidal  bump  or  indentation 
centered  about  the  body  midpoint  and  extending  over 
0.393  <_  x <_  0.607.  The  ordinates  of  the  shapes  are 
given  by 


The  next  set  of  results  displayed  in  figure  6 
are  for  a similar  shape  but  with  the  maximum  dia- 
meter located  at  x - 0.70.  The  ordinates  for  this 
shape  are  given  by 


with  n - 6.03  and  A given  by  equation  (16) . Again 
the  Cp^- input  theoretical  predictions  are  in  very 
good  agreement  with  the  data.  The  free  air  results 
indicate,  however,  as  may  have  been  anticipated 
from  the  comparisons  given  in  figure  4 and  5,  that 


bunp  where  boundary  layer-shock  wave  Interaction 
effects  may  have  occurred.  The  free  air  results 
indicate  that  tunnel  interference  exists  for  this 
shape  at  all  the  Mach  numbers  shown,  with  the 
strongest  interference  again  at  Moo  - 1.00  and  1.025. 


where  (x^,X2>  “ (0.393,  0.607) 


The  results  indicated  in  figure  7 are  for  the 
bumpy  body.  In  the  SLOR  calculation  of  these  re- 
sults, as  well  as  the  indented  body  results  to 
follow,  the  same  r-grid  used  for  the  smooth  body 
results  given  in  figures  4-6  was  employed.  The 
x-grid  was  refined  by  clustering  40  equally- 
spaced  points  between  0.3  <.  x <.  0.7,  and  then  using 
30  geometrically-expanded  points  with  grid  ratio 
1.2:1  to  cover  each  of  the  regions  -2.0  <.  x <.  0.3 
and  0. 7 <.  x < 3.0. 


The  corresponding  indented  body  results  dis- 
played in  figure  8 on  the  other  hand  show  a very 


Fig.  8 Experimental  (§)  and  theoretical  ( ) 

pressure  distributions  for  an  indented  parabolic 
arc  body  with  D/i»l/14  at  several  transonic 

Mach  numbers:  Cp^fr/D-4)  input, 

free  air 


Fig.  7 Experimental16  (§)  and  theoretical  (“ 
pressure  distributions  for  a bumpy  parabolic- 
arc  body  with  D/t-l/li  at  several  transonic 

Mach  numbers;  Cp(r/I>-4)  input, 

free  air 


small  tunnel  effect  at  every  free  stream  Mach 
number  other  than  M„  » 1.00.  The  body  indentation 
apparently  significantly  reduces  the  size  and 
strength  of  the  supersonic  pockets  for  M»  < 1,  and 
similarly,  the  subsonic  pockets  at  M„  > 1.  In 
fact,  at  K,  ■ 1.025  and  1.10  the  flow  at  r/D  = 4 is 
purely  supersonic.  We  note  that,  for  all  of  the 
cases  shown,  the  theoretical  results  for  both  free 
air  and  C_  -input  indicate  a higher  pressure  peak 
at  the  body  midpoint  than  the  data.  This  is  un- 
doubtedly due  to  a thickened  boundary  layer  in  the 
indented  region  caused  by  the  decelerating  flow. 
This  effectively  results  in  a shallower  and  milder 
indentation  and  consequent  smaller  pressure  peak. 


With  regard  to  the  results  indicated  in  figure 
7,  the  experimental  flow  field  pressures  at  the 
radial  location  r/D  - 4 — the  outermost  station  at 
which  data  were  taken — display  much  steeper 
gradients  than  any  of  the  smooth  body  data  shown 
in  figures  4-6.  In  addition,  the  data  at  H,  ■ 

1.025  and  1.10  indicate  multiple  regions  of  rapid 
compression.  Comparison  with  the  theoretical 
Cpm-input  results  displays  very  good  agreement  at 
all  the  Mach  n\xnbers.  Some  quite  minor  dis- 
crepancies are  evident  in  the  region  behind  the 
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Effect  of  Control  Surface  Location 

The  final  results  shown  in  figure  9 display 
the  effect  on  the  predicted  surface  and  near  flow 
field  pressures  of  varying  the  radial  location  of 
the  control  surface  at  which  the  experimentally- 


measured  boundary  condition  is  imposed.  These 
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Fig.  9 Effect  of  control  surface  location 
on  near  flow  field  and  body  surface 
theoretical  pressure  distributions 
at  M„=1-00  for  the  <x/-£>RmaX=0-3 
body  shape 

results  are  for  the  (x)^  = 30%  body  at  M„  = 1.00, 

for  which  comparisons  werf  previously  shown  in 
figure  5.  This  example  was  purposely  selected  for 
its  strong  tunnel  interference  characteristics  as 
well  as  the  fact  that  experimental  flow  field  pres- 
sures were  obtained  for  this  flow  out  to  the  radial 

location  r/D  = 8.  The  dashed  ( ) and  solid  ( ) 

lines  in  figure  9 correspond  as  before  to  theoreti- 
cal results  for  free  air  and  C-  -input  at  r/D  = 4. 
The  dash-dot  (-•-•)  and  dotted  (•••)  lines  corres- 
pond to  theoretical  results  for  CPm'  input  at  r/D  =2 
and  8.  For  the  results  with  C^- input  at  r/D  = 2, 
the  predicted  pressures  on  the  body  surface  and  at 
r/D  * 1 are  essentially  indistinguishable  from  the 
data  everywhere,  and  in  particular  in  the  vicinity 
of  the  shock  wave.  The  Cp^-input  at  r/D  = 4 
results  are  also  quite  satisfactory  and  differ  but 
slightly  from  the  r/D  * 2 results  only  near  the 
shock.  For  the  Cpm~input  at  r/D  * 8,  a consistent 
discrepancy  in  the  shock  positions  for  both  the 
surface  and  flow  field  pressures  that  amounts  to 
approximately  3-4%  of  body  length  is  observed, 
however.  The  flow  survey  location  at  r/D  * 8 is 
slightly  more  than  halfway  (4/7  H/£)  from  the 
tunnel  centerline  to  the  closest  point  on  the  tun- 
nel walls.  For  these  results  it  is  likely  that 
the  non-axi symmetric  effects  of  the  square  tunnel 


configuration  are  the  cause  of  the  discrepancy.  If 
the  tunnel  cross  section  were  cylindrical,  the 
radial  location  of  the  control  surface  could 
probably  be  taken  closer  to  the  wall  without  loss 
of  accuracy  to  the  interior  calculation.  The  im- 
portant point  of  these  results,  however,  is  they 
demonstrate  that,  for  any  test  model/ tunnel  confi- 
guration, there  will  always  exist  a minimum 
distance  from  the  tunnel  walls,  closer  than  which 
it  is  inappropriate  to  obtain  outer  boundary  data. 
Within  that  minimum  distance  from  the  wall, 
factors  such  as  tunnel  flow  nonuniformities,  local 
wall  disturbances,  etc.  will  tend  to  wash  out  the 
accuracy  of  the  experimental  measurements  with  a 
subsequent  detrimental  effect  on  the  theoretical 
calculations. 


Potential  Utility 

The  ultimate  utility  of  any  wall  interference 
assessment  procedure  is  in  the  evaluation  and  re- 
interpretation of  data  obtained  in  a wind  tunnel 
to  corresponding  free  air  conditions.  Toward  that 
end  there  are  several  interesting  options  in  which 
the  present  procedure  can  be  used.  For  the  results 
reported,  we  have  employed  the  assessment  procedure 
primarily  for  the  observation  of  two  separate 
effects  present  in  the  data,  viz:  (1)  nonpotential 
effects  of  viscosity  and  rotation,  and  (2)  wind 
tunnel  interference.  Nonpotential  effects  were 
evaluated  by  noting  the  difference  between  the  data 
and  the  theoretical  prediction  with  Cpm-input, 
while  the  difference  between  the  theoretical  free 
air  calculation  and  that  with  C^-input  provided  a 
measure  of  the  tunnel  interference  present.  While 
these  two  calculations  in  themselves  do  not  rein- 
terpret the  data  to  free  air  conditions,  they  serve 
to  identify  the  range  of  test  conditions  where  tun- 
nel interference  is  minimum  for  a specific  test 
article.  Such  information  is  of  major  importance 
for  existing  conventional  transonic  tunnels  where 
wall  control  is  limited  or  impossible,  and  relates 
to  "the  principle  of  minimizing  interference 
rather  than  corrections"1®. 

For  the  actual  correction  or  reinterpretation 
of  the  test  data  to  free  air  conditions,  Kemp^  has 
suggested  an  iterative  procedure  which  effectively 
corrects  the  tunnel  data  through  corrections  to 
the  oncoming  uniform  flow.  The  first  step  of  the 
method  uses  the  above  predictive  method  and  consists 
of  solving  the  potential  flow  problem  using  the 
experimentally-measured  boundary  condition.  This 
calculation  yields  the  theoretical  body  surface 
velocity  distribution  (<f>s)m.  Next,  the  correspond- 
ing free  air  potential  solution  is  obtained  which 
provides  the  surface  velocity  ((j>g) pA«  The  differ- 
ence criteria  | (4>s)m  - ( 4>s ) fa  I < e is  then  examined 
along  the  entire  surface,  and  the  data  classified 
as  negligible,  correctable,  or  uncorrectable  (a 
judgement).  If  correctable,  the  free  air  problem 
is  recomputed  in  an  iterative  sequence  which,  for 
the  axisymmetric  flows  considered  here,  consists  of 
changing  the  oncoming  Mach  number.  When  the  dif- 
ference criteria  is  satisfied,  the  tunnel  data  are 
interpreted  as  free  air  data  at  the  new  Mach  number. 
For  the  three-dimensional  case,  corrections  in 
oncoming  flow  angle  would  enter  the  iterative 
process  as  well. 
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